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A MULTILEVEL SYSTEM OP ALGORITHMS FOR DETECTING 
AND ISOLATING SIGNALS IN A BACKGROUND OP NOISE 

L. S. Gurin, K. A. Tsoy 


ABSTRACT. In a number of cases (in parti- 
cular when conducting a controlled space experi- 
ment), the problem of detecting and isolating 
signals in a background of noise should be solved 
in several stages. To start with, the information 
is processed with the help of less laborious algo- 
rithms, and then on the basis of such processing, 
a part of the information is subjected to further 
processing with the help of more precise algo- 
rithms. In this paper, such a system of algo- 
rithms is studied, a comparative evaluation of a 
series of lower level algorithms is given, and 
the corresponding algorithms of higher level are 
characterized. 


1. Statement of the Problem 

An experiment in outer space, in particular a controlled experi- 
ment, is a complex system with many levels of hierarchies from many 
aspects. As is well-known, three types of hierarchies are distin- 
guished [1]: strata (levels of description and abstraction), layers 
(levels of complexity of the given solutions), and echelons (organi- 
zational levels). Let us consider one stratum — the statistical 
processing of information. It consists of several layers and eche- 
lons; moreover, the levels are not independent with regard to these 
aspects. We are only interested in the layers. 

*Numbers in the margin indicate pagination in the foreign text. 
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Assume that In the spacecraft, results from the measurement of 
a random process are obtained which are an additive mixture of signal 
and noise . Ji These results are processed onboard by means of a specific 
algorithm, and then the processed materials are transmitted for sub- 
sequent processing on the ground by means of other algorithms. The 
purpose of the processing is to detect and isolate the signal with 
maximum accuracy and reliability. Here it is necessary to take into 
account the limitations on the high-speed response of the electronic 
computers onboard and on the ground and the carrying capacity of the 
telemetric channels. The development of a multilevel system of algo- 
rithms for the solution of this problem Is required. 

The establishment of such a system of algorithms is subdivided 
into several stages. 

1) Various models of the signals and noise are considered. To 
work out the methodology, the simplest model is chosen. 

2) Various algorithms for working out the lower level are 

considered. Their characteristics (reliability of detection, the 
time for realization on an electronic computer) In the assumed models 
of the signals and the noise are studied. It turns out that even in 
the simplest cases such an Investigation can be only partially com- 
pleted analytically. Therefore, the universal method of statistical 
modeling is applied (even though it is very laborious). /4 

The characteristics obtained depend on certain selected para- 
meters of the tuning (threshold values). These parameters are 
selected during the solution of the system problem, i.e., at one of 
the later stages; therefore, at the first stages it is desirable for 
the characteristics to achieve independence from the parameters. 

This increases the volume of the calculations investigated, as well 
as the volume of fixed results. 

3) All this is repeated for the algorithms of higher level — 
only more precisely and more laboriously. 
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4) A concrete physical problem is considered. Based on the 
processing of the first results of the measurements, conclusions are 
drawn concerning the adequacy of the models of the signal and the 
noise which were considered. If necessary, stages 1) - 3) are 
repeated for models which are in better agreement with the physical 
problem. 

5) On the basis of all of the preceding, a concrete system of 
algorithms for the given physical problem is developed, including a 
determination of the tuning parameters, etc. 


In the present paper, we discuss the first results obtained 
during the completion of stages 1), 2), and in part 3). 

2. Models of the Signal and the Noise 

In connection with the discrete character of the measurements 
and the finite length of the realization, the following form of 
representation is convenient. Given the values: 

Ui =f(? i )+£ i =s i+£ i . CD 

on the segment /«= i t . . . , A/ , where Is the signal, and 

is the noise. For various models of the signal and the noise, it is 
convenient to introduce provisional notation. The model A±. /f(x) 
is a triangular splash, i.e., 


«S f - =0 when | 1 - i 0 \ ; 

Si =c(i - h' „ hen 


( 2 ) 


are independent, normally-distributed random quan- 

t ities , where : 


+ a i X + Q 2 » 

a? = (> 0 x* + 6iX + b z . 


(3) 
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Let us consider the following particular cases of the model A-^ 

2 _2 


- O' O- -O' i k are known, i is unknown but _ 

11 t «• o p 

e terminate . rV 

.2 0® *°° 

A.«“"Ci=Cl X +Q + & , the rest similarly. Here 

12 t O 1 2 


n a a are determinate but unknown. 

2 _ 2 


2 *6 / 

A — /) . =rO ; CT =G: K is known, C and i are unknown, but 


M5 

determinate. 


A^ — the same as in A^j only k is also unknown. A-j^, A.j^, 

A l4 — the same as A ii» A 13» A l4 s respectively, but the 

unknowns (i ; i and C: k, i and C) have given laws of dis- 
0 0 o 

tribution and unknown parameters. 

The set of particular cases of the model A^ may be increased. 
,The general characteristic feature is the uncorr elatedness of the 

p 

noise. For the present we shall consider a to be known 

The model A q. In contrast to the model A^, the noise is a 

normally distributed random vector with mathematical expectation 

A. =s and covariant matrix: 

H = l&ijljur - ■ (4) 

Here, a 1 can be determined in accordance with (3). 

In the particular case (stationary random process), we obtain 
the model A^s 


o/-°i CiJ ><*(/-*) . 


(5) 


For the present we shall assume that the function (5) is known. 
The presence of a prime or a second subscript in the designation of 
the model A 2 (e.g., A’^) has the same meaning as in the case of the 


model A 


1* 
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3. Algorithms of the Lower Level 


A) The algorithm of moving summation. Let us introduce here 
the average: 


for the models 1 


i*k 


yi =-sjkiZ, y, 


( 6 ) 


j*i- k 


AAA 

i{ * x i2 > M >11 


. . as k assumes its known value. 
The question as to the models A-^, A-^, ••• we shall not consider at 
present for the given algorithm. It is clear that y . is determined 

for i A • 


To determine the acceptance of the solution concerning the 
presence of the signal or its absence, one of the following decision 
rules can be applied. 


Decision rule A^ . Let 


y-m ax = -max y i . 


(7) 


Let us introduce two decision functions (for the models A 1 and A 2 ): 


^2 = Vtnajr"*J^f|=^r » 


( 8 ) 


where , i.e., F t >F z . , and the cases f *^0 p F x* O are impossible 


The decision rule has the form: 


F. 4 o / 
•; /•; > o / 

f l >o,f 2 <o 


there is no signal 

there is a signal 

the situation is 
indeterminate 


(9) 
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It is possible to write analytic expressions for the charac- 
teristic of this algorithm. Henceforth we shall consider groups of 

n l <n points, i.e., in formula (7) the limits of variation of i 
will be still narrower. For the probability of a false alarm, we have: 


/7 


" 1 tfspoj 


2 \ 

y x P y 2 2oy * 

where the region ft in the space t^ y . . . T is 

system of inequalities: 4 

2 */ * c a + CT ; ^ = i,2 - 


( 10 ) 


the 


(ID 


The formulas for the probability of the admission of the signal and 
for the situation of indeterminacy are obtained by algorithms. Since 
the computation of integrals of high multiplicity reduces, one way or 
another, to statistical modeling, it is more convenient to apply 
statistical modeling to the direct calculation of the characteristics 
of interest to us. The results of such calculations will be dis- 
cussed below. 


Decision rule A ^. In place of (8), we introduce two functions 
of the number of the point J: 

T ; . (li 

Let us designate two further whole numbers n 1 and n 2 - The 
decision rule has the following form: 

- if the number of points of the group for which Q 

is less than or equal to n^, there is no signal; 
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- if the number v 2 of points of the group for which 
is greater than or equal to n 2 , there is a signal; 


/8 


- in all the remaining cases, the situation is indeterminate. 


For this rule it is difficult to obtain analytically the proba- 
bility of a false alarm or the admission of a signal. However, we 
shall present some estimate, or more precisely, the mathematical 
expectation of the numbers and v 2 for various groups of points 

(containing or not containing the signal). Let y denote a number 
which equals 1 or 2 . Then in the case of the model where there 

is no signal, the probability of the event f ls: 




y-2 

e~“ dt - i-OCCp) . 


(13) 


In connection with the fact that the theorem about the mathe- 
matical expectation of the sum is also valid for the independent random 
variable, we have: 




(14) 


If there is a signal and it is (has a maximum) at the point numbered 
i , then we have for y 


o 


Vi. 


Pi 




. • itk+l 2 • k(2 k+i) ’ i ’ 

; (2-k+i-j)(2k*a-j) C 


k(*k+ i) ; 

; y>aA. 


(15) 


Next , 


p (i ± i ) = I - cp (c y. ) 


( 16 ) 
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Let n i denote the set of value ±j for which the point numbered % ±j 
lies in the group of points under consideration, and let be the 
number of points of fL . Then for the model A.^ in case there is a 
signal we have: 




( 17 ) 


Decision rule A 


3 * 


This rule differs from rule A 2 in that the 


numbers and X2 replace and v 2 , respectively, where x^ is the 
maximum length of the series of events (])>' O for the group of 


points N^. It is natural that for rule A^ the values c u and n^ must 
be selected differently than for rule A^. 


For decision rule A^ also, it is possible to obtain analytic 
expressions for some estimate of the characteristics (for the charac- 
teristics themselves, this is difficult). We shall present some of 
them. 


In order to estimate the mathematical expectation of a number of 
series (in the present case we also include a series consisting of a 
single point), note that each series begins with a jump in the corre- 


sponding level, i.e., for example. 


Let us find the probability of a jump at the given point j in 
case there is no signal. When the value C is given, we have for 

r* 

the desired probability: 



where, as usual. 


n 


®(t)=7F=sl 


( 19 ) 


For the mathematical expectation of the number of series in the 
group of N 1 points, we obtain the value . Since the total number 


of points above the given level in the absence of a signal is fur- 
nished by the formula (14), then according to a familiar theorem 
of the theory of restoration, for the provisional mathematical 
expectation m of the length of a series (provided the series occurs) 
we have: 


m 




( 20 ) 


Thus, when C.. =2.C> calculation gives P: - C C2C6 , i . e . , 
r J 


m =.*1*095: Hence the provisional probability of obtaining a series 

of length greater than unity is very small. This is confirmed by 
experimental results. Note, however, that formula (20) is imprecise. 
The reason for the imprecision of formula (20) lies in the fact that 
the random variables which appear in the application of the theorem 
from the theory of restoration do not satisfy the condition of 
independence . 


The probability of obtaining a series of length x> beginning at 
the point numbered j, is given by the following formula: 


2&+x* 2 


pj u)dt { j 


( 21 ) 


V*i 


where the -region Q is defined by the inequalities: 

A 


W2 k 


2 t-c > C H a \i 2k+ i K r 2. , x * i : 


V = 

V 


( 22 ) 


r t v <Cp cri'zk*! , v c =i M v.-x + a. 

V-Xc 


As regards formula (21), we may' repeat the same reasoning which was 
stated above regarding formula (10). 
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Using the decision rules introduced, it is possible to draw the 
following conclusion relative to a data array N: if there is a signal 
and a situation of indeterminacy arises in, say, even one of the 
groups of — information is transmitted to a higher level. So far 
we have considered the model A.^. If it is necessary to evaluate 
the quantity C (in the model A^), then, depending on which rule is 
applied, we proceed in the appropriate manner. For example, for rule 
we may assume that: 


c - ik+L a 

° " {■ y max 


OWGINAC PAGE r$23) 
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In the case of the model A^» before applying the algorithm 
described, it is necessary to eliminate the trend (3) by the usual 
methods, after which we obtain the model without trend, e.g., A-^. 
Let us call such a procedure the algorithm of moving summation for a 
model with a trend . By contrast, we shall call the following algo- 
rithm the more precise algorithm of moving summation for a model 
with a trend. 


a) Let us apply the usual algorithm of moving summation for a 

model with a trend: i.e., as a preliminary measure we exclude the 

trend . 

b) If we have obtained the solution "no signal” or "indeter- 
minate", the calculation ends. 

c) If we have obtained the solution "there is a signal" — we 
estimate the parameters of the signal and subtract the signal from the 
initial process with the trend and repeat the calculation. The new 
solution is definitive. 

B) The algorithm "forward-backward". This algorithm was worked 
out and used to solve another problem concerning the subject, and it 
is discussed in [2]. It consists of the following. For each value 
of i from '£4 7 s i~k> let us draw two straight lines of 

regression, forward and backward, through the A'+l points: 


/II 
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(" 2 ) (" 2 ) /( 12 ) / • • \ • . • i 

- " 6-y ) ; ,/=%••• 2+ ^- 


v 


(24) 


Vj - a t 

J \ / 

As a statistic, on the basis of which the decision rule is constructed, 
let us take: 






/ . 


(25) 


For the dispersion o' A of the random quantity £±('i), we have [cf. 
[ 2 ]): 

.a 
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A ~fc(k+l)(k+z) 


( 26 ) 


We shall construct the decision rule by analogy with the algorithm A) , 


. w 

only we shall replace y^ by A (i ) and b y ex P ression (26). 


Decision rule . Let 


A 


m =maxA(i) . 


(27) 


Let us introduce two decision functions ( F-M): 


i max 


^ li O' / fcffc +i y 


(28) 


and apply the decision rule (9). 


Decision rule . The same as A 2 , except that in place of (12). 
we write: 


F (j)-A(j) C H ai 


(29) 


Decision rule B .,. Analogous to A^, except that the function 
(29) is used. \ 
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The advantage of the algorithm B is that in the case of the 
model A 12 , the preliminary elimination of a trend is not required. 

C) An optimal filter. It is easy to see that the algorithms 
discussed above are linear filters. If k is known exactly, then 
it is possible to use an optimal filter, giving the maximum signal to 1 
noise ratio at the output. It Is easy to show that it has the form: 

i+fi-i 


(30) 


where 


■ - Vi I Wh 
J ~ M 1 TT7 ' 


Here in formulas (8) and (12) we must write Q 


of 


Cm <T 

Formula (15) is replaced by the following: 


,a{ 


l 




3 k 


(31) 


instead 


6k>-4k-2kj+$(j 2 -^ ) 


\j 


f^(ik-j)(Zk-j-i)(zk-j+i) ; 
O : 


j=k..~ y 2k-u | 
j>2k-2. 

C k 


(32) 


For the algorithm A) the signal to noise ratio is ^ , and for 

the algorithm C) it is Q^^fSS-±^. s i. e . , we obtain an increase of 
/( 2 l)(2k+ 1 ) ' 

I'' — — ' times > which equals 1.21 times when k = 6. This 

Insignificant increase does not compensate for the substantial 
increase in the labor of calculating with C) in comparison to A). 
Moreover, for C) it is mandatory to know the parameter k. Therefore 
we did not consider the algorithm C) in more detail. 


/13 
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Note that for the algorithm B), the signal to noise ratio 
equals : 


ORIGINM- 

q£ EQQ* 9V*** hJk 


c_ Jdi+msin: 

a V Qk 



When k = 6, this is 4.-5? less than for the algorithm A). However, 
we consider this algorithm, since when k increases, it improves- -^in 
comparison to the algorithm A) (when lc - 2, it is already 1.25 times 
better). Moreover, as was pointed out above, it does not require 
the exclusion of a trend. The final concludions will be presented 
below. Finally, note that some increase in the signal to noise ratio 
for the smoothing algorithm can be obtained without an increase in 
the labor involved, due to the proper choice of the width of the 
smoothing interval. In formula (6) we shall take k^ in place of k 
(assuming it is known ) ; it also replaces k in formulas (8) - (12), 
etc. In formula (15), when J = 0, we obtain: 


t . 

' y.-~ 


t 0 2^Y**i 




Here the signal to noise ratio equals: 


(33) 


/14 


o 


V 2^,+i 


(3*0 


Let us find the maximum of expression (3*0, assuming that k^ is 
continuous and ranging over the interval (1, k) . Thus, we obtain: 


lc 


l.opt 


7 - 3* 


(35) 


In our example when k = 6, we must take k^ = This improves 
the signal to noise Indicator by 10?. 
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4. Upper Level Algorithms 


Various realizations of methods for detecting and separating 
signals from the noise (with precision with respect to optimization, 
the number of iterations, etc.) based on the study of probability 
functions, which are familiar from the literature (cf., 3-g., [3, 4 ] ) 
may be taken as upper level algorithms. For the model A-^, we have 
the following expressions for the probability functions / in the 
case of a null hypothesis (there is no signal), and •/ ^ in the case 
a signal is present: 




. . . a 


exp\ 


X> 


N 

■2 

< = 1 


L V 1 v Ui-Q-i f 


of 




( 36 ) 


1 { 2 x ) A ' /2 6 


1 


a 


exp\ 


N 


■ N 

42 


Qi)*\ 

— . 


i = i 


or: 

t 


)* 


( 37 ) 


Note that these formulas are not applicable to the model A^ . In the 
case A| and equally for the case A 2 , it is possible to obtain other 

t LI . ... 


expressions for and 
is to be changed somewhat. 


in which event the following discussion 


First of all, by whatever means the problem of detection is 

solved. It is necessary to estimate the parameters, i.e.. In essence, 

to solve the problem of isolating the signal. This has led to a large 

set of algorithms. Let us consider the method of maximum plausibility 

only, which reduces for constant a to the method of least squares. 

The maximum plausibility [ corresponds to the maximum of the 

i 

quantity: 



Depending on the model assumed, the quantity S is a function 
of the following parameters: 



- in the model 

- in the model 

- in the model 

- in the model 




A 

A 
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S = 5 (efo) ; 

S =S(i 0 ,a o ,a lf a e ); 
S = S(a,i c )i 
S - S (k, c, i 0 ) : 


ORIGINAL PAGE l§ 
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and also possible is the variant A-^ y » where 

S =S(k,Oj 0 ,a o ,a i ,a e ). 


1 

0 E P°° 


To minimize. S it is possible to apply arbitrary algorithms. The 
choice of the best of them depends on the number of unknown para- 
meters: for one parameter Cease A-^) the "golden section" method may 
be used; for two parameters, the Gauss-Seidel method may be used, etc. 
We shall make only two comments. 


In the first place, the problem of minimizing S may turn out to 
be multi-extremal, in which case one of the methods of global optimi- 
zation which is described in [5] may be employed. 


Secondly, in connection with the fact that estimating the para- 
meters is subordinate to the general problem of detection, the 
requisite precision of the estimation algorithm depends on the results 
of comparing the hypotheses. Here we have our own hierarchy of algo- 
rithms: after having estimated the unknown parameters in the first 
approximation, we solve the detection problem with the help of two /1 6 

decision functions F^ and F, as this was done at the lower level, and 
only in the event of indeterminacy do we proceed to a more precise 
estimate of the parameters. When performing the last-mentioned esti- 
mate it was possible to have used only one decision rule (excluding 
the possibility of indeterminacy), provided we did not have in view a 
further type of hierarchy (e.g., with respect to the number of 
realizations to be used). 

Let us proceed to the solution of the detection problem, assuming 
that the parameters are known. As is well knownj^ in the case of a 
single decision function, this problem is solved as follows: let us 
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find the plausibility ratio: 



( 39 ) 


and let us assign a certain threshold c. A solution concerning the 
presence of a signal is assumed if: 


t>, c . 


(Ho) 


In the contrary case, we shall assume that there is no signal. As 
regards the choice of the threshold c, this depends on the solution 
of a system of problems as a whole: the criteria, the limitations, and 
the available information. 


In the case of two decision functions (at the lower or inter- 
mediate stages of the hierarchy) let us assign two thresholds c^ and 
° 2 (c „><:,) , and let us define: 


F~£-e L ; 

Fi = t- c, . 



(41) 


Everything that follows corresponds to what was discussed above. 


For systematic calculations, it is necessary to know the proba- 
bility characteristics of the algorithms obtained. In the simplest 
cases they can be written out in explicit form in terms of the 
quantity l. Taking into account, however, the fact that the true 
values of these characteristics still depend on the accuracy of the 
preliminary parameter estimates, we shall not write down these 
formulas. The method of statistical modeling is successful in the 
case of various models and concrete algorithms in obtaining the 
characteristics, taking into account all the decisive circumstances. 

5. The Results of Statistical Modeling 

Statistical modeling was carried out on the model A 12 with the 
number of points N = 100. All together, there were 100 realizations 


2 

with a constant value for a and a random trend (3), the coefficients 
a 0 , a 1 and a 2 in the expression 

Q^-a 0 i l +a t i*a 2 

being chosen uniformly distributed on the intervals: 

a- (-0.0005; 0.0005); * Oj- (-0.05; 0,05); (^-(O^). 

■ I 

The signal was formed in accordance with (2) for £ = 50 9 k~Q > 
the coefficient C assumed three values, depending on the version of 
the calculation: C=3C ; C = Cf and C-o^CS . 

Next, when using the moving summation algorithm, the trend was 
excluded in advance with the aid of the least squares method. The 
trend was not excluded when the "forward-backward" method was used. 

Then portions of length N 1 = 25 were processed, using various 
decision rules. Here when there was a signal each of the versions 
represented 100 realizations, and when there was no signal — 200 
realizations (since in the general interval N - 100, a central por- 
tion with a signal, and two extremes without signal, were selected). 
On each portion in each realization, the values of y max (correspond- 
ing to A „ for the "forward-backward" algorithm), the number of 
max 

rejections and the maximum length of the series of the criteria 
y Cor, correspondingly. A) for the levels C equal to the following 
set of values: 2 ; 2.5 3.0 G* \ 3,5 O* and A , were fixed, 

where by a’ is meant a- for moving summation, and d. for the 
"forward-backward" algorithm. 

In Table 1 - 14 below are given the estimates of the proba- 
bilities: 

P&max^^y) » P ? P v ) * P (Xp-^^z) 

for various versions. Versions which have no interest from the point 
of view of obtaining conclusions are not presented. 
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Since the smaller values of C turned out to be more effective, 

_ M ^ 

the versions C =0,5^ CT and 0 = 50’ were calculated for y when 

=0 and C^=i,5 0* . The corresponding results are presented 
in Tables 15 - 18. 

The following experiment was performed regarding the upper 
level algorithm. In the model with a trend which was previously 
described, the trend was eliminated. Next, by the method of maximum 
plausibility i Q and C were estimated and the detection problem was 
solved according to the criterion of the plausibility ratio (40). The 
distribution of the plausibility ratio £ for the cases when the signal 
was present or absent with respect to the 100 realizations (as well 
as for the signals 0.33 cr, a and 3cr ) was outstanding for the choice 
of the threshold C. 

The results obtained are presented in Tables 19 - 20. 

The distribution h'c-501 (the true value of i is 50) is pre- 
sented in Table 19 when there is a signal for the versions C = a /19 

and C = 3cr. 

In Table 20, the distribution of £ is given in the absence of a 
signal and for three versions when it is present. 

Besides the results presented in the tables, the required 
computing time was estimated. An analysis of the results obtained 
and some preliminary conclusions are presented in section 6. 

6. Analysis of the Results Obtained and Conclusions 

Before analyzing the results obtained, it is necessary to 
comment on accuracy. As is well known, when estimating some proba- 
bility p concerning an empirical frequency which was obtained as 
the result of N trials, 5% of the confidence interval equals 

± 1 , 961 /®^^. ( 42 ) 
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In our case, when there is no signal, i.e., when N = 200, for 
probabilities of a false alarm which are nearly zero, the estimate 
turns out quite crude ; thus, for example, for p = 0.10 we obtain an 
error of + 0.025. However, account must be taken of the fact that 
when two estimates obtained by different algorithms are compared, 
the situation turns out more favorable because the trials, on the 
basis of which these estimates are obtained, are dependent, since 
all the algorithms are subjected to exactly the same realizations. 

The basis for this position is in [6] (and in earlier works of Yu. G. 
Polyak) . 

Further, if the results obtained by statistical modeling are 
compared with certain theoretical results presented in section 3, 
it must be kept in mind that the former are obtained on a model with 
a trend, and the latter — on a model without a trend. Therefore, 
the degree of their coincidence makes it possible to judge the 
quality and the influence of the trend exclusion. 

Finally, it must not be forgotten that the results obtained 
pertain to a case of uncorrelated noise and a concrete value of the 
number k. Here we are limited by the preliminary analysis, making 
it possible to obtain some conclusions with regard to the investi- 
gation of a concrete system. 

The difficulty of the preliminary analysis consists in the fact 
that the comparison must take a number of criteria into account: 
the time required to realize the algorithms on an electronic com- 
puter, the probability of a flase alarm ol., and the probability of 
the admission of a signal o^. However, in individual cases, the 
realization times for versions being compared are practically iden- 
tical. In the first place, realization occurs during the comparison 
of various decision rules and. various threshold values for one and 
the same lower level algorithm. Moreover, it occurs with sufficient 
accuracy in practice even for different algorithms of lower level, 
since the operation of trend exclusion turns out not to effect cal- 
culation time significantly. Therefore, the preliminary analysis 
of lower level algorithms can be performed without taking the time 
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'O' 

criterion into account. As regards the criteria and a 2 > f° r 

the preliminary analysis, one of them could be fixed and compared with 

the other. Unfortunately, this approach is impractical in view of 

the fact that usually the values of and a 2 do not agree for the 

versions to be compared (cf. the tables). This difficulty may be 

avoided as follows. Suppose for one of the versions we have the 

values and ct 2 , and for another, the values and a£. Let us 

consider an hypothetical system such that for the first version m 

independent detection channels are used, and the signal is assumed 

to be detected if it is detected in at least one of the.. channels; 

for the second version, let the corresponding number of channels be /21 

m’ . Then the probability of a false alarm is ma^ for the first 

system, and rn'ct'^ for the second. Let us select m and m' from the 

condition: 

T71Q& , -371 ex' ( 43 ) 

i i. 

Now it is possible to compare the systems with respect to the general 
probability of the admission of a signal: the first system is 

preferable if: 


(SXQ*rAL PAGE IS 
gl KIR QUALITY 

or 

. i 

of < cx'5 (44) 

Thus, for example, for the decision rule A 2 , let us compare 
the threshold levels C^= C)y (Table 15) and = 2(5“ (Table 2) 

when a weak signal G = 0^5bCS' is detected. In the first case (for 

1^ — 7 ) we obtain Q (l L - 0*0$ , and in the second case (for = i) 
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C3C i =:0,05. The corresponding values of a 2 are 1 = 0.26 = 0.74 
and 1 - 0.14 = 0.86. Comparison according to formula (44) shows 
that the first version is to be preferred. 

As the result of a similar analysis, the following conclusions 
are obtained. 

1. The best of the lower level algorithms is the algorithm A) 
with decision rule A 2 and C^=Q^- . 

Here the choice of the number iy is determined by the solution 
of the network problem. 

2. The previous conclusion may be changed in favor of algorithm 
C) in case: 

a) the noise is correlated; 

b) of smaller values of the number k; 

c) a trend has been eliminated in advance. 

In each of these cases additional investigation is required, 
analogous to that discussed above. 

3. According to all the indicators, except for the calculation 
time, the upper level algorithm is better than the best of the lower 
level algorithms. In terms of calculation time, it is only half as 
good. 


4. Except for the preliminary conclusions mentioned, when a 
physical problem is to be solved, the results presented in Tables 
1 - 20 may be used. 


ORIGIN AK RAGE fit 
OF POOR QUALITY) 
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TABLE 1. DISTRIBUTION 

PC iw> £ a 3 ■* 


h 

T No signal 

f 

• * 

t 

j C - 0,33© 

j 

! C«S 

t 

i 

t C 
i . 

-0,2 

1,00 

1,00 

I , Q 0 

x , oo . 

-0,1 

0,99 

1,00 

1,00 

1,00- 

0,0 

0,98 

1,00 

1,00 

1,00 

I 0,1 

* 0,91 

0,96 

1,00 

1,00 

0,2 

0,71 

0,82 

0,97 

1,00 

0,3 

0,45 

0,59 

0,88 

1,00 

0,4 

0,23 

0,36 

0,70 

1,00. 

0,5 

0.10 , 

0,18 

0,47 

1,00 

0,6 

0,006 

0,09 

0,26 

1,00 

0,7 

0,0 

0,0 

0,12 

0,99 

0,8 

- 

- 

0,05 

0,91 

0,9 

*i - - 

- 

0,0 

0,78 

1.0 

- 

- 

- 

0,63 

’ 1. 1 

- 

- • 

- 

0,51 

1.2 

- 



0,24 

1.3 

- 

- 

• - 

0,14 

1.4 



- 

0,07 

1.5 

- 


- • 

0,0 


translator' s note. Commas in numbers represent decimal points. 



TABLE 2. 


DISTRIBUTION 


P 20 




/2H 


A 1 

“7 T 

^ No signal 

C - 0 , 

336 1 C-6 
1 

j C -56 

'o / 

1,00 

1,00 

1,00 

1,00 

I 

0,05 

0,14 

0,33 

1,00 

2 

0,0 

’ 0,08 

0,26 

0,99 

3 

- 

0,05 

0,20 

. 0,99 

4 


0,02 

0,14 

0,99 

5 

- 

0,01 

0,11 

0,99 

6 

- 

0,01 

0,08 

0,99 

7 

- 

0,01 

0,06 

" 0,94 

8 

- 

0,01 

0,05 

0,92 

Q 

♦V 


v 0,01 

0,04 

0,81 

10 

- 

, 0,01 

0,02 

0,68 

II 

- 

0,0 

0,0 

0 , 41 . 

12 


- 

- 

0,24 

'13 

- 

- 

- 

0,11 

14 

- ■ 

• 

- 

0,01 

15 

- 

- 

- 

0,01 

16 

- 

- 

- 

• 0,0 


^Translator ' s note. Commas in numbers represent decimal points. 
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*Translat or ’ s note. Commas in numbers represent decimal points 
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TABLE 5. DISTRIBUTION P (%>*») FOR c * = 2,5 S'r * 

I 

OkiGINAi. f AGE IS 

PB BOOK QUAUTTC 


h : 

t No signal 

1 

! C = 0 . 33 S 
; 

! 0 = 6 " ! 
i i 


. 0 

1,00 

1,00 

1,00 

1,00 

i 

0,0 

• 

0,01 

0,12 

0,99 

* 

2 

- 

0,0 

0,09 

0,97 

• 3 

- 

- 

0,05 

0,93 

4 

\( 

- 

0,02 

0,88 

5 

- 

- 

0,02 

0,81 

6 

- 

- 

0,02 

0 , 75 ' 

7 

/ 

t 

0,01 

0,73 

8 

- 

- 

0,01 

0,65 

9 

- 

- 

0,01 

0,51 

10 

- 

- 

0,0 

0,36 

II 

- 

- 

- 

0,17 

12 

- 

- 

- 

0,07 

13 

/ 

- / 

- 

- 

0,01 

14 

- 

- 

0,0 


"Translator’s note. Commas in numbers represent decimal points. 
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TABLE 6. DISTRIBUTION P(\‘p=s-£U) EORCm^SSu* 


£ ' ! 


f I No signal, 


! 04=0,33.5* ! C = 6 ! C =5^ 


1 ,00 


1 .00 


1,00 

0,04 

0,01 

0,01 


1,00 

0,90 

0,85 

0,78 

0,70 

0,64 

0,58 

0,48 

0,29 

0,16 

0,07 

0,04 

0,01 

0,01 


^Translator r s note. Commas in numbers represent UwT 


TABLE 


7. DISTRIBUTION P ^ &y, ) I f 0R. C^- S^Cj* * 


U, . 

1 No Signal 

t 

i 

! C= 0,336 

T 

j 

! C =G ' 

! 

t 

• -- • J » 

! C =$5 

\ 

i 

0. 

1,00 * 

1,00 

1,00 

1,00 

I 

0,0 

0,0 

0,04 

0,90 

2 • 

- 


0,01 

0,83 

3 


- 

0,01 

0,7o 

4 


! 

i 

0,0 

0,67 

5 

«■ 

- 

- 

0,60 

6 

* 

- 

V 

0,54 

7 


1 

I ^ 

- 

0,42 

8 


- 

- 

0,26 

9 



- 

0,15 

10 


- 

- 

0,07 

II 


- 

- 

0,04 

12 


- 

- 

0,01 

13 

- 

- 

- 

0,01 

14 

- 

- 

- 

0,0 • 


^Translators note. 


Commas in numbers .represent decimal points. 









TABLE 9. DISTRIBUTION P(\^> h) FOR Cy ~ 2 * /3I 



^Translator ’ s note. Commas In numbers represent decimal points. 
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TABLE 11. DISTRIBUTION P( v V ^ j for (V-2.5G 


h 

? No signal! 

f 

! C = O.SSS' 

t 

! C 

t 

0 

1,00 

1,00 

I 

X 

0,28 

0,30 

0 

2 

0,08 

0,14 • 

0 

3 

0,03 

... # ' 

0,03 

0 

4 

0,05 

0,01 

0 

5 

0,005 

0,0 

0 

6 

0,0 

•*> 

• 

7 

• 


• 


! C *6 L 0*55* 

t » 


TABLE 12. DISTRIBUTION P( ty, ) FOR C 


0 

? . No signal 

\ 

! C = 0,33 6 

» 0 = 6 “ 

! C =3 S' 

i 

• 

0 

' 1,00 

1,00 

1,00 

1,00 

I 

0,28 

0,29 

0,47 

0,97 

2 

0,05 

0 r 06 

0,16 

0,87 

3 

0,01 

0^01 

0,06 

0,73 

4 

" 0,0 

AO 

0,01 

| 0,50 

5 

• .•* • " 

/ - 

0,0 

0,16 

6 

■T : . 

/ ‘ T 

** 

0.0 


•Translator's note. Commas in numbers represent decimal points 








TABLE 13. DISTRIBUTION FOR A,' 3 $ 


o 

j No signal 

! C a 0,33 & 

! C=s S'* ! 


tv 

.! 

f 

t 

T ~t 

. . . i . * 

« X 

o 

1,00 

*'*••■** II u ■’ * "*■' ■- * 

2,00 

1,00 

1,00 . 

■' I 

0,095 

0,n 

0,20 

0,86 . 

2 

0,015 

o.ex 

0,87 

d,?l 

3 

0,05 

0,01 

0,02 

©,4? 

4 ‘ 

• 0,0 

0,0 

0,0 

®,25 

5 




0,@2 

6 


. » 

/ 


©,0 


TABLE 14, DISTRIBUTION POL* W FOR C 


«. 1 

No signal 

X 0=0,33 6^ 
/! 

/ i 

! C =<S ! 

i i 

•• • 

C=3 (S' 

0 

1,00 

1,00 

1,00 

1,00 

I 

0,095 / 

' 0,11 

0,20 

0,86 

2 

0,01 

0,01 

0,07 

0,68 

3 

0,005 / 

0,01 

0,01 

0,46 

4 

0,0 

/ 

. 0,0 

0,0 

0,24 

5 

i 

/ 

«K» 

- 

0,02 

6 . 

- / 

/ 

/ 

f 

— 

— 

0 6 () 


"Translator's note. Commas in numbers represent decimal points 
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TABLE 16. 
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DISTRIBUTION P(V^v) FOR V *.56g 


ti 

? No signal 

| 

f 

t 00,336 

| i 

! C= $ 
! 

j 

! C » 

t • . 

i 

. 0 

1,00 

1,00 

1,00 

1,00 * 

I 

0,52 

. 0,60 

0,86 

1,00 

2 

0,36 

#■ 

0,48 

0,?9 

1,00 

3 

0,24 

• 0,41 

0,70 

1,00 

4 

0,16 

0,36 

C .59 

1, 00 

5 

0,12 

0,27 

0.54 

1,00 

: 6 

0,10 

0,22 

0,50 

I, GO 

! 7 

0.07 

0,17 

0,44 

1,00 

| 8 

0.04 

0,10 

. 0,36 

1,00 | 

9 

0.02 

0,08 

0,29 

1. 00 

: 10 

0,005 

0,06 

0,22 

0,99 

! n 

0,005 

0,05 

0,18 

0,97 

i 12 

0,005 

0,04 

0,16 

0,94 

13 

0,005 

0,03 

0,14 

0,72 

14 

0,0 

0,0 

0,07 

0,39 

15 

- 

- 

0,04 

0,23 

16 

- 

— 

0,02 

0,10 

17 

- 

- 

0,01 

0,02 

18 ' 

- 

- 

0,0 

0,0 


^Translator ’ s note. Commas in numbers represent decimal points. 


34 










fcttGINAt 
•®® *°°R to 


18. DISTRIBUTION POL* 2 j; ) for Cm « 1,5 6. j 


" X ! 

No signalj 

J 0 = 0 , 336 " 

t 

| C 

T 

! Os 3 

i : . 

• 

/ 0 

1,00 

I.oo 

1,68 

1,00 

I 

0,19 

0,25 

0,67 

1,00 

2 

0,12 . 

0,20 

0,51 

1,00 

3 . * ' 

* 0,06 

0,12 

0,37 

1,00 

4 

0,02 

0,08 

a , 35 

0 , S 9 

5 

0,02 

0,30 

0,30 

0,99 

6 

0,01 

0,04 

0,26 

0,99 

7 

0,005 

0,03 

0,17 

©, S 9 

8 

o 

« 

o 

0,02 

0,12 

0,99 

9 . 


0,02 

0,09 

0,97 

10 

- 

0,01 

0,08 

0,91 

ii _ ... 

, - - ' 

0,01 

0,07 

0,83 

12 

- 

0,0 

0,03 

0,58 

13 . 

- 

- 

0,02 

0,32 

14 


- 

0,01 

0,17 

15 


- 

0,0 

0,03 

16 

* *v 

«* 

"■ ' 

0,0 


*Translat or ' s note. Commas in numbers represent decimal points. 
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PC l* L) 


/39 


* 


C «= 0,336 

! c = 6 

t 

? 

* Cass' 

( 

• 

f 

1,00 

1,00 

1,00 

0,94 

0,99 

1,00 

0,76 

0,97 

1,00 

0,61 

0,93 

1,00 

0,45 

0.81 

1,00 

0,36 

0,66 

I, GO 

0,26 

0,53 

1,00 

0,17 

0,43 

1,08 

0,13 

0,36 

1,00 

0,07 

0,23 

0,93 

0,04 

0,14 

0,97 

0 t 0I 

0,06 

0,79 

0,0 

0,04 

0,59 


0,01 

0,40 


0,0 

0,22 



0,15 



0.0 





